Activated aging dynamics and negative fluctuation-dissipation ratios 
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In glassy materials aging proceeds at large times via thermal activation. We show that this can lead 
to negative dynamical response functions and novel and well-defined violations of the fluctuation- 
dissipation theorem, in particular, negative fluctuation-dissipation ratios. Our analysis is based 
on detailed theoretical and numerical results for the activated aging regime of simple kinetically 
constrained models. The results are relevant to a variety of physical situations such as aging in 
glass- formers, thermally activated domain growth and granular compaction. 
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Glassy materials display increasingly slow dynamics 
when approaching their amorphous state. At the glass 
transition, where relaxation times exceed experimentally 
accessible timescales, they change from equilibrated flu- 
ids to non-equilibrium amorphous solids. In the glassy 
phase, physical properties are not stationary and the sys- 
tem ages m . A full understanding of the non-equilibrium 
glassy state remains a central theoretical challenge. 

An important step forward was the mean-field de- 
scription of aging dynamics for both structural and spin 
glasses In this context, thermal equilibrium is never 
reached and aging proceeds by downhill motion in an 
increasingly flat free energy landscape • Two-time cor- 
relation and response functions depend explicitly on both 
their arguments, and while the fluctuation-dissipation 
theorem (FDT) does not hold it can be generalized using 
the concept of a fluctuation-dissipation ratio (FDR) . This 
led in turn to the idea of effective temperatures jjL and 
a possible thermodynamic interpretation of aging 0, |^ ■ 

However, in many systems of physical interest, such 
as liquids quenched below the glass transition or do- 
main growth in disordered magnets, the dynamics is 
not of mean-fleld type, displaying both activated pro- 
cesses and spatial heterogeneity 0, 13 • While some ex- 
periments and simulations nonetheless seem to de- 
tect a mean-field aging^egime, theoretical studies have 
found ill-defined FDRs [3 , non-monotonic response func- 
tions 10, in [Hill, observable dependence [Hill, non- 
trivial FDRs without thermodynamic transitions and 
a subtle interplay between growing dynamical correlation 
lengthscales and FDT violations 0, experiments 
have also detected anomalously large FDT violations as- 
sociated with intermittent dynamics It is thus an 
important task to delineate when the mean-field concept 
of an FDR-related effective temperature remains viable. 

Independently of this interpretation, FDRs have addi- 
tionally been recognized as universal amplitude ratios for 
non-equilibrium critical dynamics. This makes them im- 
portant markers for distinguishing dynamic univers ality 
classes. This area has seen a recent surge of interest j20j . 



Here we study kinetically constrained models, which 
are simple models for glassy systems with heterogeneous 
dynamics . We use them to study systematically the 
impact of activated, and therefore strongly non mean- 
field, dynamics on FDRs and associated effective tem- 
peratures. In addition, the dynamics of these models 
becomes critical at low temperatures, where dynamical 
lengthscales diverge. Our work therefore also pertains to 
the study of FDRs in non-equilibrium critical dynamics. 
We show that FDT violations retain a simple structure 
with well-defined FDRs in the activated regime, elucidate 
the physical origin of negative dynamical response func- 
tions, and predict the generic existence of negative FDRs 
for observables directly coupled to activated processes. 

We focus mainly on the one-spin facilitated model of 
Fredrickson and Andersen (FA model), defined in 
terms of a binary mobility fleld, Ui € {0, 1}, on a cubic 
lattice; Ui — 1 indicates that a site is excited (or mobile). 
The system evolves through single-site dynamics obey- 
ing detailed balance with respect to the energy function 
H = Si '^2- The dynamics is subject to a kinetic con- 
straint that permits changes at site i only if at least one 
nearest neighbor of i is in its excited state. In any spatial 
dimension d, relaxation times follow an Arrhenius law at 
low temperatures |^ . 

At low temperatures, T < 1, the dynamics of the FA 
model is effectively that of diffusing excitations, which 
can coalesce and branch ^2§, . Such a problem can 
be described in the continuum limit by a dynamical field 
theory action with complex fields Lp{r,t), >p{r,t) p3.l25|. 



Cp{dt - DV"^ ~ 7)(p - 7(^V + + (1) 



with Hi — > {1 + ip)ip |25j| . In terms of the equilibrium 
concentration of excitations, c = (rij) ~ 1/(1 + e^/-^), 
the effective rates for diffusion, coalescence and branch- 
ing are, respectively, D oc c, A oc c, and 7 oc l23|. 
The aging dynamics following a quench from T = 00 to 
low T consists of two regimes. Initially, clusters of exci- 
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tations coalesce on timescales of 0(1), reaching a state 
made of isolated excitations. This process and the sub- 
sequent energy plateau are reminiscent of the descent to 
the threshold energy in mean- field models "^l. For times 
larger than 1/D, excitations diffuse via thermal activa- 
tion and decrease the energy as they coalesce. 

Following we calculated the connected two-time 
correlation to all orders in A at tree level, 



Here Cq{t,t^) is the Fourier transform of 
{ni(t)nj{t^))c = {n^{t)nj{t^)) - {ni{t)){nj{t„)), 
and ipq{t),(pq{t) are those of (p{r,t),(^{r,t). Subscripts 
q indicate wavevector, N is the system size, and 
z = Dq^tvr- We have assumed that both waiting time 
and final time t are large {':^ and set 7 = 

to get the leading contribution at low T. The function 
f{z) goes as / w 1 — 1/z for z 3> 1, and / w (1 + z)/3 
for <C 1. At g = we get the energy auto-correlation: 
Co(t,tw) - N-^H{t)H{t^))^ « n{t)t^/{3t), where 
n{t) — {ni{t)) « is the mean energy per site. 

These classical expressions should be accurate above the 
critical dimension dc where fluctuations are negligible. 
The limit 7 = corresponds to diffusion limited pair 
coalescence (DLPC) which has dc — 2 (25i | . 



Consider now a perturbation SH 



-hqAq at time 



with Aq — J2i cos(<7 • ri)ni. The action changes by |27| 
S S f f 

OlT-q Jp^s Jp 



to leading order in c. For the two-time response function 

NRq{t,t^) = TS{Aqit))/6hq{t^) = -T _q{t) {t^)) 

we then find at tree level 

Rqit, iw) « (At2)-ie-^9'(*-*w)(z - 1). (2) 



The energy response follows as i?o(i,iw) ~ —n{t)/t. The 
corresponding susceptibility, xo(^:iw) — /j dt'Ro{t,t') 
is always negative, and proportional to the density of 
excitations at time t, Xoit^tw) ~ ""-(0(1 ~ tw/t). 

The FDR Xq{t,t^), defined through Rq{t,U) = 

Xq{t, tvj) dt„Cq{t, tv 



if T is included in Rq, reads 



Xq{t, t^) 



Z-1 



{l + d,)zf{z) 



1-1/z (z>f) 
-3 + 12z (z < 1) 



At any waiting time t„, for inverse wavevectors q^^ 
smaller than the dynamical correlation length, ^(tw) — 
\/Dt„, FDT is recovered: Xq 1. However, at small 
wavevectors, q^ <C 1, the FDR becomes negative. In the 
q ^ limit, one gets the FDR for energy fluctuations, 



Xoo = ^o(i,tw) = -3 (d>4). 



(3) 



This simple form means that on large, non-equilibrated 
lengthscales a quasi-FDT holds with T replaced by 
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FIG. 1: Normalized FD-plots for the FA model at T = 0.1 
in d = 3 with fixed t — 2 x 10^ and running € [0,t]. 
Symbols are results from simulations, full lines are the tree- 
level predictions. Wavevectors shown are q — tt/x with x — 1, 
2, 2.4, 3, 3.2, 4, 5.33, 6, 8, 12, and 16 (top to bottom). System 
size was iV = 32^. Avera ees are over 4 x 10* initial conditions. 



T/Xoo- Contrary to pure ferromagnets at criticality, 
however, Xoo is negative. This feature is not predicted by 
earlier mean-field studies as it is a direct consequence of 
thermal activation: if temperature is perturbed upwards 
during aging, the dynamics is accelerated; the energy de- 
cay is then faster, giving a negative energy response. 

The fluctuation-dissipation (FD) plots in Fig. 1 show 
that the tree-level calculations compare very well with 
numerical simulations of the FA model in d = 3. We have 
also confirmed the diffusive scaling with z = Dq^t^. Sim- 
ulations were performed using a continuous time Monte 
Carlo algorithm. We measure Cq{t,tvi) as the auto- 
correlation of the observable Aq, on a cubic lattice with 
periodic boundary conditions, N — and a linear size 
L ^ \/Dtv/. The susceptibility Xq{t,t^) is obtained by 
direct generalization of the "no-field" method of jl^ to 
continuous time. We show data using the prescription 
of Rcf. 29], plotting Xq{t,'tw)/Cq{t,t) as a function of 
1 — Cq{t,t^)/Cq{t,t) for fixed observation time, t, and 
varying waiting time, < < t. The abscissa runs 
from (tw = t) to ~1 {tw <C t). Using as the running 
variable ensures that the slope of the plot is the FDR 
Other procedures, e.g. keeping fixed, 
would give very different results 0, 0, O, 0, 01 that 
could lead to erroneous conclusions. 

In dimensions d < dc 'we need to take fluctuations into 
account. For d = 1 exact scaling results can be derived 
for the FA model by considering a DLPC process with 
diffusion rate D = c/2. Because the long-time behavior 
of DLPC dynamics is diffusion controlled 12^ we are free 
to choose the coalescence rate as A = I? without affecting 
scaling results in the long-time regime {t, ^ 1/c) }3(|. 
We focus on the low temperature dynamics (c — > 0) and 
times t<Cl/70cl/c^ where branching can be neglected; 
the system is then still far from equilibrium as it ages. 

To analyze our DLPC process we use the stan- 
dard quantum mechanical formalism [Mj: probabil- 
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FIG. 2: Normalized energy FD-plots for the FA model in 
dimensions d = 1 to 4, for fixed times t and running waiting 
times t„ G [0,t]. Symbols show simulation data, for T = 0.2 
to 0.06 and t — 10"* to 10^. Full curves show a straight line 
with slope Xoo ~ —3 and the exact limit plot in d = 1 with 
asymptotic slopes shown as dashed lines. 



FIG. 3: Normalized FD-plot for the local observable Ui in 
the d = 1 FA model at T = 0.1 and t = 2 x 10*, 5 x 10*, 
lO'^ and 10*^ (bottom to top). The dashed diagonal is FDT. 
Inset: final part of the t = 10^ data, enlarged to show the 
non-monotonic response and Xoo ~ —3.307 (full line). 



ities are mapped to states \P{t)) — 
where \n) = |ni, . . . , tiat), so that the master equa- 
tion reads dt\P{t)) = Wc\P{t)), with Wc the DLPC 
master operator. Correlations are then {A{t)A{tvr)) = 
|/l|^gWcAt^gWct„|p(o)), where (1| = and At = 

t — tw DLPC dynamics is closely related to diffu- 
sion limited pair-annihilation (DLPA) [sj by a similar- 
ity transformation B: Wc — BWaB~^ where Wa is 
the DLPA master operator with diffusion rate D = c/2 
and annihilation rate 2A = c. Introducing empty and 
parity interval observables, Ei ~ nlLii(l ~ "-fe) and 
Pi = rifcJiill - 2nfe) respectively, where i = (11,^2), 
one has {l\Eie^'^^^ = (l|Pie^^^*B"\ {1\B = (1| and 
{l\niB — 2(l\ni [sj. If we interpret the Ui in DLPA as 
domain walls in an Ising spin system, rii = ^{1 — atai+i), 
then Pi = (7^2 , and Wa = becomes the T = mas- 
ter operator for the Glauber-lsing spin chain. Substitut- 
ing the two-spin propagator (IjcTij o'i^e'^"'^'^* derived in 
[33| and mapping back to DLPC then gives 

(l|£;,e^-^*=i/,,_,,(2cAt)(l|+ ^ Gg(cAt)(l|i?,-, 

jl<j2 

(4) 

where Gg(-) and i?„(-) are given m |33]. Eq. Q to- 
gether with the identity rii = \ — E(^i^i^i) is the key in- 
gredient for the calculation of the two-time energy cor- 
relation and response functions Co(t, iw), ^0(^7 ^w) in the 
d — 1 FA model j^. Qualitatively a picture similar 
to d > dc emerges, but with the FDR now showing a 
weak dependence on the ratio t/t^. At equal times, 
Xoit/t^ ^ 1) = -(1 + V2) = -2.414..., while for 
t/tvj — > 00 the FDR approaches 

Xoo ==37r/(16-67r) -3.307... id=l). (5) 

Fig. 2 demonstrates complete agreement between our 
scaling predictions and simulation data in all dimensions. 



In d = 1 the data fall on the curved limit FD-plot |30| 
derived from In dimensions d = 3,4 the simulations 
are compatible with a constant FDR Xoo = — 3, Eq. 
The data for d = dc = 2 suggest a slightly curved limit 
plot, possibly due to logarithmic corrections, but are still 
compatible with Xoo ~ — 3. 

Numerically, the no-field method of l28l| becomes unre- 
liable for small q. To get precise data as shown in Fig. 2 
we used instead the exact relation 

2i?o(i,tw) = (l - 2c)dtn{t) + dt„Coit,t^) + C4t,t^). 

for the energy-temperature response of a general class 
of kinetically constrained spin models |3Cll |. The quan- 
tity C*(t,tw) is defined as N-^H{t)U{t^))c, with U = 
Si ~ c) and fi € {0, 1} the kinetic constraint. 

The auto-correlation C{t,t^) = {ni{t)ni{t^))c, and as- 
sociated auto-response, i?(t,tw) = T5{ni{t)) /5hi{t^), to 
a local perturbation 5H = — hiri i , a lso have a negative 
FDR regime. Previous studies [l3,ll3 bad suggested that 
the corresponding FD-plot has an equilibrium form, even 
during aging. A more careful analysis, however, reveals 
nonequilibrium contributions (Fig. 3). Exact long-time 
predictions for C{t,t^) and R{t^t^), and the resulting 
FDR, can be derived from Q 30]. One finds that, as for 
the energy, the FDR has the aging scaling X{t/t^). It 
crosses over from quasi-cquilibrium behavior, X « 1, at 
t/tyj K 1, to X Xoo for t/tw ^ 1; notably, Xoo here is 
the same as for the energy, Eq. However, the region 
in the FD-plot that reveals the nonequilibrium behavior 
shrinks as 0(l/\/t) as t grows. This explains previous 
observations of pseudo-equilibrium and makes a numer- 
ical measurement of Xoo from the local observable rii 
very difficult. On the other hand, Fig. 2 shows that with 
the coherent counterpart, i.e. the energy, this would be 
straightforward, confirming a point made in [T5l |. Numer- 
ical simulations for d > 1 produce local FD-plots analo- 
gous to Fig. 3 but with Xoo = —3. We emphasize that 
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FIG. 4: Normalized energy FD-plots for the East model for 
(T,t) = (.15,100), (.15, 10^), (.2,50) and (.2,5000) (from left 
to right). Inset: Normalized FD-plot for local observables for 
T = 0.15 and f = 10^, 10^ and lO'^ (from bottom to top). 



the non-monotonicity observed is not an artefact of using 
t as the curve parameter, as in e.g. |^.IT^. 

We have shown that some important components of 
the mean-field picture survive in systems with activated 
dynamics, in particular the concepts of time sectors (ini- 
tial relaxation versus activated aging on long timescales) 
with associated well-defined FDRs. However, activation 
effects make these FDRs negative. This effect has not pre- 
viously been observed in non-equilibrium critical dynam- 
ics, and calls into question whether effective temperature 



descriptions are possible for activated dynamics. Our 
arguments show that negative non-equilibrium responses 
should occur generically for observables whose relaxation 
couples to activation effects; activation need not be ther- 
mal, but can be via macroscopic driving, e.g. tapping of 
granular materials 0,|3^. To illustrate how quantitative 
effects may vary, we consider briefly the East model 
of fragile glasses, leaving all details for a separate report. 
The behavior is richer due to the hierarchical nature of 
the relaxation, which leads to plateaux in the energy de- 
cay. Properly normalized energy FD-plots nevertheless 
have a simple structure, with three regimes (Fig. 4): for 
given equilibrium FDT is obeyed at small time dif- 
ferences t ^ t^, indicating quasi-equilibration within a 
plateau; a regime with a single negative FDR follows, 
coming from the activated relaxation process preced- 
ing the plateau at t; finally, the plot becomes horizon- 
tal, corresponding to all previous relaxation stages which 
decorrelate the system but do not contribute to the en- 
ergy response. Interestingly, in the FD-plots for auto- 
correlation and auto-response (Fig. 4, inset), each relax- 
ation stage of the hierarchy is associated with a well- 
defined effective temperature, a structure reminiscent of 
that found in mean-field spin glasses j^. 
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